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A stack of annular detectors made of high purity germanium and a magnet spectrograph were 
used to measure pp — > dir + differential cross sections at a beam momentum of 850 MeV/c over a 
large angular range. A total cross section of a — 0.2301 ± 0.0036 (stat) ± 0.0230 (syst) mb and an 
anisotropy A2/A0 = 0.856 ± 0.016 were deduced. These values follow fits to low energy data. From 
the present A2 value it is found that the pionic p-wave amplitude ao is larger than assumed so far. 

PACS numbers: 



I. INTRODUCTION 

Although the reaction pp — ► dn + has been inten- 
sively studied, some open problems remain. These 
are mainly connected to the importance of differ- 
ent partial wave amplitudes. For small pion ener- 
gies only three complex amplitudes are believed to 
contribute to the physical observables: ao, ai and 
d2- Here the index denotes the angular momentum 
in the proton-proton channel. Close to threshold, 
the cross section is given by only one amplitude 
a\ corresponding to an s-wave in the final chan- 
nel. For slightly larger energies, the p-wave am- 
plitudes ao and ai show up. Because \a-i\ is much 
larger than |ao|, the latter amplitude is often ig- 
nored @. The d-wave amplitudes are even smaller. 
In the simplest experiments, three observables can 
be measured, the total cross section, the anisotropy 
and the asymmetry of the reaction products. The 
usual method to deduce the three complex ampli- 
tudes is then to make use of the Watson theorem 
Q. However, this procedure often yields negative 
values for |a | ||. Because of its smallness, the 
magnitude of ao cannot be extracted from the to- 
tal cross section. It enters the other two quantities 
via interference with other amplitudes. Especially 
in the anisotropy, there is an interference between 
the two p- waves. Thus the measurement of this 
quantity might yield an estimate of ao. 

The energy dependence of the anisotropy, which 
will be defined later, is a puzzle. Data published 
before 1996 like those of Ref. [Ef lie on but of- 



ten below the predictions of the SP96 phase shift 
analysis carried out by the Virginia grou p |5|| . Near 
threshold data, appeared first in 1996 || Qand are 
in agreement with the SP96 solution. The newest 
data are from a measurement of the inverse reac- 
. They are even below the Ritchie et al. 
, thus indicating the negligible effect of ao 
even in the interference term (see section III). This 




is in contrast to SAID calculations as well as to ex- 
trapolation of the new near threshold results. Fur- 
thermore, the new data of Ref. § show a much 
larger scattering than the threshold data, although 
the latter are from two different experiments. It is 
worth mentioning that in the measurement of the 
anisotropy, the luminosity cancels out. Large fluc- 
tuations in this quantity may thus point to system- 
atic errors other than target thickness and beam 
intensity. We, therefore, investigate the situation 
by a measurement in the range of these data not 
so close to threshold. 



In the present experiment, the differential cross 
section for the pp — > dn + reaction was measured 
at a beam momentum of 850 MeV/c, which cor- 
responds to a pion centre of mass momentum di- 
vided by the pion rest mass of rj — 0.51. This 
value is within the range of the Pasyuk et al. data 
Q and will thus allow to proof the importance of 
pp-wsive interference in the anisotropy as already 
mentioned. 



II. EXPERIMENTAL PROCEDURE 



The measurement of the reaction pp — ► dir + was 
performed using a detector that combines large 
momentum and geometrical acceptance for heavy 
recoiling reaction products. A proton beam with a 
momentum of 850 MeV/c was extracted from the 
COSY accelerator and focused onto a target cell 
containing liquid hydrogen. It had a diameter of 
6 mm and a thickness of 6.4 ± 0.3 mm [|l0| with 
windows of 1.5 fiin Mylar. The excellent ratio of 
hydrogen to heavier nuclei in the window material 
reduced empty target events to a negligible level. 
The beam spot had dimensions of less than 1 mm 
and divergences of better than 4 mrad. Beam halo 
events were suppressed using a plastic scintillator 
as a veto counter with a 4mmdiameter inner hole 
in front of the target. The detector system is the 
"Germanium Wall", which is part of the GEM- 
detector at COSY in Julich |11 and the remod- 
elled magnetic spectrograph Big Karl . Here we 
give only some additional details specific for this 
experiment. 

The magnetic spectrograph was set to zero de- 
grees in the laboratory system. It was used to mea- 
sure the low energy branch of deuterons close to 
the primary beam. The germanium wall consists 
of three high purity germanium detectors with ra- 
dial symmetry with respect to the beam axis. The 
first detector (called Quirl-detector) measures the 
position and the energy loss of the penetrating par- 
ticles. The active area of the diode is divided on 
both sides into 200 segments by 200 grooves. Each 
groove is shaped as an Archimedes' spiral covering 
an angular range of 2n. They are mainly used for 
measuring the energy loss of the penetrating parti- 
cles or the total kinetic energy of stopped particles, 
respectively. The next detectors are divided into 
32 wedges to reduce the counting rate per division 
leading to a higher maximum total counting rate 
of the whole detector. 

Fig. [I] shows the response of the germanium 
wall for reaction particles. Clearly distinguished 
are two bands corresponding to protons from the 
pp — > ppir and deuterons from the pp — > dir + re- 
actions. The quantities measured with the germa- 
nium wall are energy, emission vertex and particle 
type. They were converted to a four-momentum 
vector. These measurements and the knowledge of 
the four momenta in the initial state yield the miss- 
ing mass of the unobserved pion by applying con- 
servation of momentum, energy, charge and baryon 
number. 

The germanium detectors have holes in their 
centres. The primary beam passes through these 
holes and is then led via an exit in a side yoke of the 
first dipole magnet of the magnetic spectrograph 
to a well shielded beam dump. Recoiling deuterons 
at emission angles inside this hole, i. e. close to 
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FIG. 1: Energy loss in the first calorimeter crystal 
Ei versus the loss in the second named E$ (A possi- 
ble crystal E2 was not mounted in this experiment). 
The visible bands are due to detected protons and 
deuterons. 



zero degree, were detected by the magnetic spec- 
trograph. Details of such measurements are given 
inRef. §. 

The reaction deuterons were selected in the off- 
line analysis by applying gates to the kinemati- 
cal loci in Fig. [l. The efficiency of the analy- 
sis procedures were studied by Monte Carlo calcu- 
lations |Q. Finally, the data were corrected for 
reduced efficiencies due to nuclear absorption in 
the detector material |D|. The correction factors 
vary from 17% to 22% for the present energy inter- 
val. The present setup has full acceptance for the 
recoiling deuterons, except for a small area close 
to the beam exit hole. For particles not stopped 
in the germanium wall, only energy loss measure- 
ments are made. Because the angular distribution 
of the reaction products must be symmetric to 90 
degree in the centre of mass system, it is sufficient 
to measure only one half of the distribution. The 
deuterons emitted backward in the centre of mass 
system are all stopped in the germanium wall. We 
restrict ourselves to this part of the distribution in 



order to not introduce ambiguities when unfolding 
the backbending part in Fig. |l|. 

Each of the measurements, i. e. emission angle 
= (9, <j)) and energy of the recoiling deuteron, is 
sufficient to extract the angular distribution. Be- 
cause of this overconstraint, one can clearly dis- 
tinguish between reaction and background events. 
This is shown in Fig. 0. The missing mass of 




FIG. 2: Missing mass spectrum. The experimental 
result is shown by solid line histogram, the assumed 
background by dashed line histogram. 

the unobserved 7r + is shown as deduced in the off 
line analysis. The distribution has a resolution 
of 5.9 MeV/c 2 . The logarithmic scale was cho- 
sen to make the small background visible. The 
background, shown by a dashed curve in Figure |^, 
was measured with an empty target, normalized to 
integrated beam intensity and then subtracted. 

Finally, the counts were converted to cross sec- 
tions by normalizing to the target thickness and 
number of beam protons. The latter were mea- 
sured with calibrated luminosity monitor counters. 
In the calibration procedure the direct beam in- 
tensity and the corresponding number of scattered 



particles are compared. The former was measureci 
with the trigger hodoscope in the focal plane of the 
magnet spectrograph. This number is of course 
much larger and lead to dead time in the ho- 
doscope. The beam intensity was then reduced 
by debunching the beam between the ion source 
and the cyclotron injector. For sufficiently small 
beam intensity the relation between monitors and 
hodoscope is linear. The counting rate in the mon- 
itors was in the the production runs small enough 
to have dead time effects on a negligible level. This 
procedure yields a systematic uncertainty of 5% 
for the beam intensity. This yields together with 
the target thickness uncertainty of also 5% a max- 
imal total systematical error of 10% for the overall 
normalization of the differential cross section when 
both errors are added linearly. 

III. RESULTS 

The measured angular distribution in the centre 
of mass system are shown in Fig. ||. The point 
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FIG. 3: The measured angular distribution is shown 
by the full symbols. Shown are the statistical errors 
only. The overall normalization of the cross section 
has total systematical error of 10%. The solid curve 
is a Legendre polynomial fit up to second order, the 
dotted curve the SP96 solution 

near cos(#) = — I was measured with the magnetic 
spectrograph. It has larger statistical error than 
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all other points which were measured with the ger- 
manium wall. This is due to its small acceptance. 
The data show a very strong anisotropy when com- 
pared to those close to the threshold || (?J. This 
is an indication of the p-wave strength and thus 
of the importance of the intermediate A-p system. 
The angular distribution was fitted by a series of 
Legendre polynomials P 2 l[cos(9)) 
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A a + A 2 P 2 [cos(e)} + AiP4cos(6)]. 



(1) 



The fitted parameters A 2 l are connected with am- 
plitudes a,i 



A 



«2 



C(\d\ 2 )} (2) 



where C(|d| 2 ) represents all d-wave contributions, 
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with D{\d\ 2 ) denoting other d-wave contributions 
and E{dd) terms with interferences between two 
d- waves. A4 is given by 



M = - 
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The amplitude a\ is the s-wave in the pion- 
deuteron channel, ao and a 2 denote p-waves and 
03 to et6 denote d-waves. In Eqs. [| to || the no- 
tation of Mandl and Regge ji| is used. The total 
cross section can near threshold be written as 



a rj ± U\r\ 



(5) 



The first term corresponds to the s-wave amplitude 
ai and the second to the two p-wave amplitudes a n 
and a 2 . 

Fits of Eq. [l] up to second and to fourth or- 
der were performed. Although the x 2 -values of 
both fits are similar, the large uncertainties in A4 
and in A 2 favour the second-order fit. Applica- 
tion of F-statistics indicates the same conclusion. 
The results are independent whether the point at 
cos(#) = —1 is excluded or not. The negative sign 
of A4 is in agreement with a theoretical prediction 
|14| as well as with phase-shift analyses [0, |) . On 
the other hand, corresponding fits to data in this 



energy range [[| yield positive values. From phase^ 
shift analysis it is known that all terms in d-waves 
are small. Only di(ae) is of some size. If we neglect 
all other terms one gets 



A 4 = 
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18(a 3 <) + 10V5(a 4 a*)j } .(6) 



The first term is positive and the second nega- 
tive. This yields a cancelation making the re- 
sulting value even smaller in agreement with the 
present result. For the sake of simplicity we rely 
on the second order fit. This choice may intro- 
duce a systematical error. However, the change in 
the quantity of interest, i. e. the ratio A 2 /Aq is 
small; both results agree with each other within 
error bars. The second order fit is also shown in 
Fig. H as solid curve. Since A is the total cross 

TABLE I: Fitted Legendre polynomial coefficients 
(Eq. 0) in fib 



parameter 


second order 


fourth order 


A 


230.1 ± 1.2 


228.5 ± 1.5 


A 2 


197.0 ±3.1 


190.8 ±5.1 


A 4 




-6.02 ±3.9 



section we have from the second order fit a value 
a = 230.1 ± 3.6 {stat) ± 23.0 (sys) fib with the sys- 
tematical uncertainty discussed above. The effi- 
ciency correction is angle dependent, varying from 
17% to 22% in the presently measured range. If 
a 10% precision of this correction is assumed, it 
adds 1.5% to the uncertainty in the total cross sec- 
tion only. This uncertainty is added in quadrature 
to the fit error and is contained in the statistical 
uncertainty. We compare the present angular dis- 
tribution with a prediction from the SAID phase 
shift analysis ||, shown as dotted curve in Fig. ||. 
The prediction is always « 10% below the data. 
This is just the sum of all systematic errors of the 
present experiment. Whether this is the reason 
for the disagreement or not will be discussed be- 
low. From the result in the Table an anisotropy 
A 2 /A = 0.856±0.014 (fit) ±0.0045 (effic. corr.) 
is obtained. This is close to the SAID prediction of 
0.89 indicating the good reproduction of the shape 
of differential cross sections. 

This ratio is compared in Fig. ^ with previous 
results |l5|, [l6| [lTj] . The point joins the data in the 
upper part of the band formed by the experiments 
and is in agreement with the extrapolation of the 
low energy results assuming a dependence A 2 oc r\ . 
However, it is 15% larger than the Pasyuk et al. 
result. This seems to be a large discrepancy since 
only the shape of the angular distributions are 
compared. Also shown in Fig. || is the prediction 
from fit E in Ref. fl8|| under the assumption of no 
interference at all: A 2 /Aq — a\rf / (aor] + a\rf ) 
(compare Eq. 0). 
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FIG. 4: Deduced ratio A2/A0 as function of the di- 
mensionless pion momentum. The presently deduced 
value is indicated by the arrow. The earlier data are 
indicated by different symbols. The SAID calculation 
is shown as solid curve, the calculation with the as- 
sumption of no interference as dashed curve. 

It should be mentioned that the quality of fits for 
the total cross sections for the phase shift analysis 
and fit E from Ref. JlSflis similar, except for the 
near threshold region. However, the corresponding 
data are in the fit E but not in the SAID analysis. 
We have also inspected the results of the phase 
shift analysis solution C500. It only differs very 
slightly in the range 77 < 0.35 from the solution 
SP96 and gives almost identical results for larger 
77-values. It is worth mentioning that in Eq. |^, 
according to SAID, the term ^(aoa^) is several or- 
ders of magnitude larger than all the others, except 
the I d2\ 2 term. 

IV. DISCUSSION 

An almost complete angular distribution for the 
reaction pp — > dn + at a beam momentum of 850 
MeV/c was measured. This corresponds to a pion 
centre of mass momentum 77 — p^/m^ — 0.51. The 
measurement was performed with a solid state de- 
tector with axial symmetry for larger angles and a 
magnetic spectrograph for emission angles close to 
the primary beam. The angular distribution shows 
a large anisotropy when compared to data closer 
to the threshold. This is an indication that p-wave 
emission is the dominant reaction mechanism. Be- 
cause d-wave amplitudes are much smaller than 



the s- and p-wave amplitudes one cannot extract 
their strength from the total cross section. The an- 
gular distribution is more sensitive to these waves. 
The present data give almost no evidence for a d- 
wave contribution to the angular distribution. Pre- 
cise measurements at somewhat higher energies are 
desirable to investigate the importance of d-waves 
esp. the sign of the Legendre coefficient A4 in the 
range below the resonance. The ratio A2/A0 does 
not demand d-waves. However, they contribute to 
the scattering asymmetry via interferences with p- 
waves having larger amplitudes. The weakness of 
d-wave amplitudes was also recently found in the 
measurement of spin transfer coefficients even at 
400 MeV beam energy jlj| . 

The present measurement confirms a larger ra- 
tio A2/AQ as extracted from recent measurements 
at lower beam momenta and is in good agreement 
with the phase shift analysis of the Virginia group. 
It disagrees with the results of Ref. [pi. However, 
these data show a large scatter which is surpris- 
ing, since the comparison is made on a relative 
base. It has its origin mainly in large fluctuations 
of the A2 values. The agreement between the new 
data and the SP96 solution, which does not include 
these data, is much better than the earlier phase 
shift analysis Q. The present total cross section 
a = Aq and the second Legendre coefficient A2 are 
larger than the SP96 solution but the ratio A2/A0 
is in excellent agreement with this solution. A non 
negligible value for ao, as is applied in the SP96 so- 
lution, is confirmed by the present measurement. 
One needs more observables than the presently 
measured ones to deduce the amplitudes. Since 
the present result is in agreement with the SP96 
solution we will rely on this solution. It yields for 
the ratio |ao|/|a 2 | a value of 0.099 close to thresh- 
old and 0.096 for the present beam momentum. 
This can be compared to the result obtained by 
Ref. . From cross section and polarisation mea- 
surements, which, together with the above cited 
Watson theorem, a ratio |cio |/|ct2 1 = (9±36)*10~ 3 
was extracted , a value which is compatible with 
zero. 
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